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Abstract 

We apply path integrals to study nonequilibrium work theorems in the context of Brownian 
dynamics, deriving in particular the equations of motion governing the most typical and most 
dominant trajectories. For the analytically soluble cases of a moving harmonic potential and a 
harmonic oscillator with time-dependent natural frequency, we find such trajectories, evaluate the 
work-weighted propagators, and validate Jarzynski's equality. 



* Electronic address: |daveminh@gmail.com| 
^Electronic address: adiba@mail.nih.gov 



1 



As nonequilibrium work theorems describe exact relationships between path averages and 
equihbrium thermodynamic properties , path integrals offer a natural framework for 
their analysis. Previous studies of such theorems have invoked path integrals primarily with 
the intent of proving their validity, e.g. by establishing an exact relation between the path 
integration measures of processes taking place in opposite directions j^, 4|, or by invoking 
the Feynman-Kac formula for the work- weighted density 5|, |6(]. In either case, the theorems 
can be proved without explicitly computing some of the quantities that make path integrals 
most useful, for example most likely trajectories and propagators. 

Here we complement these efforts by using path integrals to obtain trajectories and 
work-weighted propagators, and explicitly evaluating these objects for analytically tractable 
models. We derive the equations of motion governing two important trajectories (Eq. flTUl) ): 
the most typical, which has the highest probability; and the most dominant, which con- 
tributes most significantly to the Jarzynski average. Interest in these trajectories stems in 
part from the observation that when typical and dominant trajectories diverge, the Jarzynski 
free energy estimate converges slowly [7]. Path integrals also offer a convenient framework 
for the computation of work-weighted propagators. We show that these propagators are 
analytically soluble for harmonic oscillators with moving equilibrium centers or changing 
frequencies. Collectively, our results demonstrate that the utility of path integrals in the 
context of nonequilibrium work theorems goes beyond simply proving their validity. 

We begin our analysis by formulating a path integral representation of the propagator, 
which for simplicity we take to be that of a one- dimensional overdamped Brownian particle 
moving on a time- dependent potential U{x,t). We assume a constant diffusion coefficient 
D and measure energy in units of fc^T. Accordingly, the propagator p = p{x,t\xQjO) from 
the point a;o at t = to the point x at time t satisfies the Smoluchowski equation [8j, 
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As with the case where U is time-independent [8|, we can apply the substitution p = 
^-u{x,t)/2q^ where q = q{x,t\xo,0), to write the equation in a Schrodinger-like form, 

Q- (2) 



dq ^d'^q 
— = D — - + 
dt dx"^ 



The "effective potential" V is given by 



V(x,t) = D 



U"{x,t) fU'{x,t) 



(3) 



2 



with primes denoting partial derivatives with respect to x. Expressingthe propagator in the 
form of Eq. ([2]) allows us to directly apply the Feynman-Kac formula [2] to write its solution 
as a path integral, 

q= [ ' dx{t)e~^^'^i^~^~^^\ (4) 

Jxo,0 

where the density q = r|xo,0) is at time t = t, V = V{x{t),t), and analogously for U. 
(In this notation, the Wiener integral over all trajectories that start at a;o at t = and end 
at X at t = r is obtained by setting V = U = in the above equation, cf. JqI). Returning to 
the original variable p = p{x,t\xq,0) and imposing the boundary condition p{x,0\xq,0) = 
6{x — Xq), we arrive at the desired path integral representation of the propagator, 

p = e-^ [ ' dx{t) e-^o^'^*[l|-v-|W] ^ (5) 

where AU = U{x,t) — f/(xo,0). Note that our approach differs from previous treatments 
3|, |4| in that no stochastic integrals (e.g. Stratonovic or Ito) are used; through the Feynman- 
Kac formalism, only ordinary time integrals and the usual Wiener contribution appear in the 
path integral. This derivation thus avoids a common confusion associated with the definition 



of the action due to such stochastic integrals jlO|. 

To make contact with Jarzynski's equality, consider the work-weighted propagator p^ 
p^(a;, r|xo, 0), defined as 

p^ ^ e-^ rdx{t) e-/o^'^*[l§-v-i^] ^_[,(,)] 

Jxofl 



e 



^ rdx{t)e-^^<^-''^'^^\. (6) 

Jxofl 



This propagator differs from Eq. ([5]) in that each trajectory in the path integral is endowed 
with the additional weight e""''^*^*^^, where w[x(t)] = j^dt{dU/dt) is the work along the 
trajectory x{t) The path average in Jarzynski's equality, (e"*") = Z^/Zq, can be written 
in terms of Pw as 

/r ^-U{xo,0) 
dxo J dx — Pu,{x,t\xo,0), (7) 

where Zt = J dye~^^'^''^\ Thus, knowledge of the analytical form of p^j reduces the compu- 
tation of the Jarzynski average (e"*") to simple quadrature. 

A unique feature of path integrals is that they allow one to extract noteworthy trajectories 



from a given path average, such as Eq. ([7]). This is done by first identifying the action t 
dictates the trajectory probability and then optimizing the desired functional (see 



hat 
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for a similar approach in the context of time-independent potentials). In the case of Eq. ([7j), 

the path average can be written as 

r r /"^''^ p-s[xit)] 
{e-"") = dxo dx dx{t)^ e-'"!"^], (8) 

J J Jxo,0 ^0 

where the action functional is defined as 



sm] ^ i^i-ii±^ + /;rft [1^ - V - If J . (9) 

Notably, in contrast to the usual case where the end-points Xq and x are fixed, here the 
action embodies the probability of observing the end-points themselves; upon minimization, 
this gives rise to the "natural" boundary conditions below. Two trajectories of interest 
are the most typical {xT{t)) and the most dominant {xo{t)). In the present formalism, 
they correspond to the trajectories that minimize the action S alone and that minimize 
the combined action-work functional, S + w, respectively. Carrying out the functional 
optimization of S and S + w with free end-points, we obtain the Euler-Lagrange equations 



x{t) _ d 
'2D ~ ~dx 



2 dt 



(10) 



subject to the natural boundary conditions a;(0) = D{dU/dx)^^^ and x{t) = 
—D {dU/dx)^^^. The sign of dU/dt in Eq. ( fTOl) is positive for the most typical and negative 
for the most dominant trajectories, whereas the boundary conditions are the same for both 
types. 

The path integral formulation of the work- weighted propagator (Eq. ([6])) and the Euler- 
Lagrange equations for the most typical and most dominant trajectories (Eq. (fTOj) ) are the 
central objects of interest in this paper. We now demonstrate them on two analytically 
tractable problems. 

(a) Moving oscillator, U{x,t) = k{x — vtY/2. Here the spring constant A; > is time- 
independent, and the center of the harmonic potential moves at constant velocity v. We first 
compute the most typical and most dominant trajectories. The Euler-Lagrange equations 
above become 

x = {Dkf{x-vt)±Dkv, (11) 
with boundary conditions x(0) = Dkx{0) and x{t) = —Dk{x{r) —vt). The sign convention 
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FIG. 1: Illustration of Eqs. (jl2p -(jl3p. The most typical (blue, bottom two curves) and most 
dominant (red, top two) trajectories for two different velocities are plotted with respect to the 
center of the harmonic well, vt (dashed). The parameters are vt = D = k = 1. 



is the same as Eq. ( |T0|) . Solving these elementary differential equations, we find 

XTit) = - ^ (1 - e-^'') , (12) 
XB(t) = t;t + ^(l-e-^'=(--*)). (13) 

Thus, in the allotted time interval < t < r, the most typical trajectory starts out aligned 
with the center of the well and lags behind at later times, while the most dominant one 
starts ahead of the center in the beginning, but falls back at the final time r (see Fig. [1]). 

We now proceed to the computation of the work-weighted propagator. For the moving 
harmonic potential, Eq. gives 



Pn 



AU I DkT 

e 2 2 X 



(14) 



This path integral can be simplified with the change of variables z{t) = x{t) ~vt — which 
does not change the integration measure. This gives 



where zq = xq — z = x — vt ~ and G is the Gaussian path integral 

G{z,t\zo,0) = [ ' dz{t) 



4£)T" 4 ^ 



(15) 



(16) 



zo,0 



This is the propagator of a Brownian particle subject to a harmonic, time-independent 
source/sink term. Its result is well-known (cf. 9[), namely 

fc (z^+Zq) cosh{DkT)—2zzo 



G{z,t\zo,0) 



exp 



sinh(Dfcr) 



^ sinh(DA;r) 



(17) 



Finally, using this result in Eq. (fT5|) and expressing the final formula in terms of the conve- 
nient variables z and zq (see above), we obtain 



exp 
Vw = 



-Uz-zo) 



(18) 



In accordance with expectations, when f = 0, this expression reduces to the well-known for- 
mula for the propagator of a Brownian particle in a stationary harmonic potential (Ornstein- 
Uhlenbeck process, see e.g. [3|), as in this case the work is identically zero and we must 
have p = Pw Additionally, using this propagator in Eq. ([7]) and performing the two Gaussian 
integrals leads to (e^"^) = 1, which confirms the Jarzynski equality prediction. 

We leave this section with a comparison of our results with those of Mazonka and Jarzyn- 



ski These authors directly solved a Fokker-Planck equation for the joint work-position 
density f{y,w\yo), where in their study y(t) = x(t) — vt. Integrating over the work, the 
marginal density of y{t) was found to be Gaussian with a mean moving according to a sim- 
ple expression (Eq. (14a) of that reference). Incidentally, upon Boltzmann-averaging over 
the initial condition y^, this reduces to our Eq. f|T2l) . We note, however, that this coincidence 
between average position and most typical trajectory is not to be expected in general, and 
is likely to be a peculiarity of harmonic problems. Lastly, we note that Eq. ( ITSl) can in 
principle be obtained from the full expression for f{y,w\yo) in the appendix of Ref. 13] by 
integrating out w from / with weight e""". 

(b) Time-dependent spring constant, U{x,t) = k(t}x^/2, where k(t) > 0. This problem 
has been used in the numerical examples of Ref. [l5|. Perhaps unsurprisingly, the most 
typical and most dominant trajectories are found to be Xr(t) = ^^(t) = for any k{t). This 
simply refiects the parity symmetry of the problem, and indicates that other properties, 
such as the width of the instantaneous position distribution, are responsible for the observed 
convergence properties of the Jarzynski average |l^. 

A less trivial aspect of this problem is its work- weighted propagator. From Eq. (j6]), we 
obtain 

p. = e-^+f J;'^*^ X rdx{t)e-^^^''i^^'^^'], (19) 

Jxo,0 

where k = k{t). Due to the fortuitous combination k + Dk"^ in the effective harmonic po- 
tential, this propagator can be computed in closed analytical form for arbitrary k{t). To 
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proceed, we factor out the Xn, x dependence from the path integral with the usual substitu- 



tion for harmonic potentials 



16|, x{t) = x(t) + y{t). The trajectory x{t) is an extremum of 



the action in Eq. (JT9l) . satisfying the Euler-Lagrange equation 



X 



{Dk + D'k')x 



(20) 



with boundary conditions a;(0) = xq and x{t) = x, i.e. y{0) = y{T) = 0. After integration 
by parts and changing the path integration variable from x{t) to y{t), we get 



e-^+^Io'^ii'e~^\o G{t) 



(21) 



where we have defined the xq, x-independent quantity 



G{r) 



0,T 



dy{t) e 



4n 4 y 



(22) 



0,0 



This Gaussian path integral can be formally evaluated in different ways. Using the Gelfand- 
Yaglom method jo!], we find G{t) = 1/ A/47rD(/9(r), where ip{t) satisfies the same differential 
equation as x{t) (Eq. ( l20l) ). with boundary conditions v?(0) = and 0(0) = 1. 

Thus, we have reduced the evaluation of the path integral in Eq. (|T9l) to the solution of a 
single differential equation (Eq. (120|) ). There are two independent solutions to this equation, 
namely 



^i(t) 



(23) 
(24) 



The desired quantities x{t) and Lp{t) are obtained by a linear combination of ipi and ^/'2 with 
the appropriate boundary conditions. Here we omit the details of this computation and 
simply quote the final result for p^, obtained by plugging these solutions into Eq. (|2T|) : 



2Dk{T) + ^) x2 - (2Dfc(0) - Ixxo] } , (25) 



where for brevity of notation we have defined ipi = ipiir) and tp2 = (''")• Thus, for any 
k{t), the work- weighted propagator is a Gaussian function whose time-dependent coefficients 
can be found via quadrature, i.e. via Eqs. fl23|) - flM|) . As in the previous problem, one can 
check that this result reduces to the propagator for a one-dimensional Ornstein-Uhlenbeck 
process when k{t) is a constant. Furthermore, upon using this expression for in Eq. ([7]) 
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and carrying out the two Gaussian integrals, massive cancelation of terms ensues, leaving 
(e"*") = >y k{0) / k{T) , which also validates the Jarzynski equality prediction. 

Lastly, we consider work- weighted propagators and dominant/typical trajectories in the 
context of forward and reverse processes, in the sense of Crooks [2]. For a given time- 
dependent energy function U{x, t) defined between times t = and t = t, let us arbitrarily 
call the dynamics under the potential U{x, t) the forward process, and the dynamics under 
the time-reversed potential U{x, t — t) the reverse process. By definition, the forward prop- 
agator = {x,t\xo,0) is given by Eq. ([5]), while the reverse propagator p^{x,t\xq,0) 
is given by the same expression with the replacement U{x,t) U{x,t — t). Now consider 
the quantity p^ = p^{xo,t\x,0), i.e. the reverse propagator with inverted initial and final 
conditions, x and xq respectively. After the time-reversal change of variables s = t — t and 
y{s) = x{t — s) in the path integral, we obtain 

p^ = e^ r"rfy(t)e-^"^^*[^-^-^^^], (26) 

where Af/ is defined as before. We thus see that p^ differs from p^ simply in the signs of 
Af/ and dU/dt. 

Two interesting consequences are of immediate notice. First, the following relation is 
verified: 

e-c/{xo,o)pF^^^ r|xo, 0) = e'^^^' V(a:o, r\x, 0), (27) 

where the forward work- weighted propagator p^ is given by Eq. ([6]). When U is 
time-independent, this result reduces to detailed balance, i.e. e~^*^^°V(x, r|xo, 0) = 
e~^*-^-'p(xo, t|x, 0), as in this case w = and the process directions are immaterial. This 
identity can thus be seen as version of detailed balance in the context of time- dependent 
potentials, and establishes that the computation of pw is no easier or harder than that of 
p. In contrast, when p^j is averaged over the initial Boltzmann distribution (by dividing the 
above equation by Zq and integratingover xq), the resultant density has the generic closed- 



form analytical solution e~^*^^''^'*/^o P, ll5|. Further integration over x leads to Jarzynski's 
equality. 

Secondly, as per Eq. ( ITOl) and the definition of reverse process above, the time-reversal 
of the typical and dominant trajectories in the reverse process satisfy precisely the same 
Euler- Lagrange equations and boundary conditions of the dominant and typical trajectories 
in the forward process, respectively. For example, observing Fig. [1] from right to left, the 
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typical trajectories become the dominant ones, and vice-versa. Our results thus provide 
means to explicitly compute and illustrate the conclusions of Ref. 

In summary, we have offered a detailed path integral study of the Jarzynski equality. This 
path integral approach allowed us to derive exact results for work-weighted propagators, as 
well as analytical expressions for typical and dominant trajectories in model problems. As 
argued above, such propagators are generally not known in closed analytical form, and 
hence one needs to resort to models such as the present ones for additional insight. Our 
Euler- Lagrange equations of motion make it possible to investigate quantitative properties 
of typical and dominant trajectories formerly discussed at a qualitative level. A natural 
application of such equations can be found in the context of single-molecule "pulling" ex- 
periments, where it can reveal the typical and dominant rupture events inherent to these 
processes [g]. 
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